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A model recently introduced for diluted magnetic semiconductors by Berciu and Bhatt (PRL 
87, 107203 (2001)) is studied with a Monte Carlo technique, and the results are compared to 
Hartree-Fock calculations. For doping rates close to the experimentally observed metal-insulator 
transition, a picture dominated by ferromagnetic droplets formed below a T* scale emerges. The 
moments of these droplets align as the temperature is lowered below a critical value Tc<T* . Our 
Monte Carlo investigations provide critical temperatures considerably smaller than Hartree-Fock 
predictions. Disorder does not seem to enhance ferromagnetism substantially. The inhomogeneous 
droplet state should be strongly susceptible to changes in doping and external fields. 



PACS numbers: 75.50.Pp, 71.30.+h, 75.40.Mg 

Recent advances in experimental techniques have al- 
lowed for the introduction of magnetic elements into 
semiconducting hosts beyond the solubility limit. This 
has stimulated the exciting research area of diluted mag- 
netic semiconductors (DMS) with the prospect of ma- 
nipulating the charge as well as the spin degrees of free- 
dom, and its possible technological applications in such 
novel fields as "spintronics" . Prominent among these 
new compounds is the III-V system Gai-^Mn^As, where 
ferromagnetic (FM) transition temperatures as high as 
110 K were obtained.El Although some of the basic in- 
gredients of the physics of these materials such as the 
local antiferromagnetic (AF) exchange between the Mn 
spins and the charge carriers are known, the origin of 
such fairly high critical temperatures is still in question. 
The ferromagnetism itself can be understood as being 
carrier-induced, in a similar fashion as the FM state 
in the double-exchaage (DE) model for manganites at 
intermediate doping.EI However, alternative descriptions 
based on Ruderman-Kittel-Kasuya-Yosida (RKKY) in- 
teractions between the impurity spins are claimed to 
lead to transition temperatures comparable to those 
observed. Q Yet it is not obvious whether RKKY can apply 
for systems with small Fermi energies as is the case for 
compensated semiconductors. Additionally, the effects 
of disorder that accompanies chemical doping was ne- 
glected in this approach.Q Recent annealing studies have 
shown the relevance of disorder and defects in DMS.eI It 
was proposedtl, based on Hartree-Fock (HF) calculations, 
that disorder in the dopant position and hopping ampli- 
tudes has a strong effect on the magnetic properties of 
these compounds, and should be incorporated in trying 
to understand the behavior of DMS. 

The above described subtle interplay between localiz- 
ing and delocalizing tendencies in these models often can- 
not be described correctly in a mean- field picture, but 
it requires nearly exact methods such as Monte Carlo 
(MC) techniques. A well-known example are the man- 
ganites in the low-doping regime, where mean-field ap- 
proximations and exact solutions lead to different con- 



clusions about the nature of the metal-insulator transi- 
tion (MIT). The cluster formation and percolative pic- 
ture that has emerged for these systems cannot-.be cap- 
tured properly with mean- field approximations a There- 
fore, here the model introduced in Ref.6 for DMS in the 
low-doping limit describing electrons in an impurity band 
is investigated using MC methods, and the results will be 
compared to mean-field calculations. 

Although the system under study is actually hole- 
doped, it is described in terms of an electron-doped ma- 
terial. Therefore, the acceptor level is treated as a donor 
level lying approximately 100 meV below the conduction 
band. Furthermore, only the impurity band itself will be 
considered, whereas the valence (conduction) band states 
are neglected. This should be a good approximation as 
long as the Fermi energy of the "electrons" is consider- 
ably smaller than the gap between the donor state and 
the lowest conduction level. The simplest such modelEl 
takes into account the hopping between random Mn sites 
and the AF interaction between the impurity spins and 
the charge carriers. The disorder potential and Coulomb 
interactions are neglected and the model is given by 

H DM s = ^2(Uj ~ Vkj){c\ a cj a + H.c.) 

+ J ^"( C ]a^ cr c t ,/3Cj / 3)S l , (1) 

where cj is the creation operator for an electron at the 
impurity site i with spin a, and iy is the hopping ampli- 
tude between sites i,j. For two sites at distance r=\i— j\ 
an exponential form, <(r)=2(l + r/os) exp(— r/as) Ry, 
was assumedB, with as the Bohr radius associated with 
the impurity site and Ry its corresponding binding en- 
ergy E b . Si represents the (assumed classical) spin at 
the Mn (3d 5 ) site with |Si|=5/2, cr denotes the vector 
of Pauli matrices, \i the chemical potential, and Jij de- 
scribes the non-local interaction between impurity sites 
and mobile holes. Owing to the localized nature of the 
Mn spins, it can be expressed as Jy= Jexp(— 2r/as), 
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with J an AF coupling constant. J is chosen as the en- 
ergy unit and periodic boundary conditions (PBC) are 
imposed. The effects of an external field can easily be 
introduced via a Zeeman term. 

Following Ref.6, we assume the hole binding en- 
ergy £^=112.4meV (lRy) and, additionally, we have 
J=15meV. It was checked that our conclusions below do 
not depend crucially on the exact values of the exchange 
coupling and, therefore, the results should be typical for 
the systems considered. The lattice constant a_pf the 
GaAs lattice is known to be 5.65A, and aB—7.8Aa With 
these values, the various nonlocal couplings and hop- 
pings tij can easily be calculated once a random distri- 
bution of the impurity ions within the host lattice (of 
size L 3 ) has been chosen. In addition, these systems 
are presumably strongly compensated through As anti- 
site donors resulting in a relatively small effective hole 
concentration of 10% of the Mn ions.EJ 

The previous mean-field study of Hums resulted in 
some interesting conclusions about the magnetic behav- 
ior of DMS.B In particular, an increase in Tc by a factor 
of ^2 was reported when a disorder Mn distribution was 
studied, compared with a periodic distribution. This is 
a challenging counterintuitive prediction that deserves to 
be tested with techniques beyond the HF approach, even 
more so as it is known that critical temperatures typically 
are overestimated by mean-field methods, as it happens, 
e.g., in the DE modcl.El Fortunately, Hamiltonian ([[]) is 
especially suited for MC calculations, and the procedure 
used here has been extensively discussed in the study of 
the manganiteso 

The MC calculations were performed on a variety of 
lattice sizes, but the bulk of our work was done study- 
ing systems with A=80 impurity sites and a doping level 
of xwO.02, which lies in the neighborhood of the realis- 
tic MIT.I Other lattice sizes (iV=50-100), with slightly 
different values of x, yield similar results and lead to the 
same conclusions. The density p of charge carriers is fixed 
at p=0.1 (with respect to the number of Mn dopants), 
unless otherwise noted. The random placement of Mn 
ions in the GaAs host lattice creates regions of high and 
low impurity density, respectively. Some of the results 
presented below are for only one such configuration of 
Mn impurities, but it was confirmed that the results for 
other configurations are within 10% of the data shown. 
Furthermore, for the configuration selected the HF crit- 
ical temperature Tq F is the same as for larger systems 
(7V=200, 300) with similar values of x,p and J. Most in- 
teresting is the magnetization behavior, with and without 
applied magnetic field. It is analyzed by measuring the 



magnetization per site, m=M/N (M=yJ2ij S» 1 Sj), of 
the Mn spins, neglecting the contribution of the charge 
carriers. The latter are of much less importance as they 
have a spin of 1/2 only, and also a density that is only a 
fraction of that of the Mn spins. Fig.[JJ compares m for 
the MC as well as for the HF solution, with both results 
obtained for the same configuration of Mn impurities. 
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FIG. 1. Magnetization of the Mn spins vs. T for model (|l]) 
with a;«0.02 and iV=80 impurity sites. Shown are results of 
MC as well as HF calculations, for either the random system 
or for impurities in a superlattice (^periodic) configuration. 
Disorder strongly increases Tc in the mean-field description, 
but much less so in the MC approach. The strong influence 
of a magnetic field at low temperatures is demonstrated in 
the inset. This is in contrast to HF, where the effect of h is 
remarkable at all T's (open circles, ft=0.08). 

In this particular case-T^^ ^0.85 (similar to earlier 
results on larger systemsl) for the impurity configura- 
tion considered, whereas the MC solution suggests an 
enhancement of FM correlations at ~0.3(±0.1). Simi- 
lar values and comparable deviations are found for other 
configurations and system sizes as well (FigJ|(a)), which 
supports the notion that mean-field treatments gener- 
ally overestimate Tc's. At this same temperature, M 
(scaled by |S|) starts to deviate from the totally disor- 
dered case M—l as well and also a pronounced change 
in dM/dT is observed (Fig.§)(b), signalling the develop- 
ment of FM correlations. Furthermore, the concept of 
a disorder-enhanced Tc does not appear in agreement 
with MC. Fig.|l| points out the influence of disorder as it 
also includes HF and MC results for the periodic case. 
From the magnetization data as well as the magnetic 
fluctuation, xoc(M 2 ) — (M) 2 , a critical behavior around 
X* er ~0.25±0.05 can be deduced for the periodic system, 
which is slightly below the critical region for the disor- 
dered case. This appears in agreement with previous HF 
resultsO, but the reduction in Tc is much less pronounced. 
It will be argued below, however, that the actual transi- 
tion temperature is lower than T^ er . A thorough analysis 
reveals that the disordered system breaks apart into oc- 
cupied and almost empty domains (in both approaches) 
and the main contribution to m is entirely due to a sin- 
gle, large cluster, with average particle density (n) cl >p. 
This leads to an increased temperature where FM corre- 



2 



lations start building up compared to the periodic case 
since the effective local magnetic field is enhanced in this 
region. Single sites or very small clusters are virtually 
unoccupied and do not order at any finite temperature. 
For this reason the magnetization is never fully saturated 
in the MC approach as in the HF description, but reaches 
a plateau at about f 5% of the possible value for a sample 
with the given Mn doping rate. If x is selected consider- 
ably lower, it would be very unlikely for such a cluster to 
exist and there should be no finite saturation maffneti- 
zation, as it is observed for samples with x<0.005.S The 
conclusions above are also supported by analyzing the 
spin correlations between different clusters A, B (FigJ^). 
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FIG. 2. (a) shows m vs. T for two system sizes, and in the 
case of the iV=80 system, the magnetization function for a dif- 
ferent impurity configuration. In all cases, 0.02<a;<0. 025. (b) 
M (scaled by |S|) for AT=50, 80 and dM/dT (triangles), (c) 
Schematic representation of the development of FM clusters 
and subsequent alignment as T is lowered. Ordered domains 
form at T* , and align at the critical point Tc- 

This reveals that_each cluster turns FM as the tem- 
perature is reducedtj, but the overall magnetization will 
remain low as the magnetization vectors for different do- 
mains point in different directions, whereas in HF they 
are required to be aligned along a chosen S^-axis. In 
other words, the idea that derealization of charge carri- 
ers leads to a-,spin alignment of the Mn ions across the 
whole samplefl, appears to be an artifact of the approxi- 
mation employed. It also leads to the unphysical obser- 
vation of totally aligned impurity spins in the extremely 
diluted limit a;<0.001 at very low temperatures. But it is 
interesting to note from Fig.|^ that the spin correlations 
between sites of a given cluster and all other sites as- 
sume a finite value at temperatures below ~0.08 (fOK). 
In addition, finite correlations develop between A, B at 
T^0.2. The exact value of this temperature likely de- 
pends on the relative position between two specific do- 
mains, but is in good agreement with the value Tc=0.08 
quoted above, where all sites were taken into account. 



Additionally, M/|S| is proportional to the number of 
spins below this temperature as it should be for a FM 
ordered state (FigJ|(b)). In any case, these values are 
lower than T c er and they suggest two different temper- 
ature scales in the problem: T* (~ 0.3), where isolated 
magnetic domains are formed, and a smaller true critical 
temperature Tc (~ O.f) where long-range order is estab- 
lished. Below T* , spins in the vicinity of a given cluster 
start to align, i.e. the magnetic domain grows in size as 
the gain in kinetic energy outweighs the loss in entropy. 
Further derealization of the charge carriers then causes 
the alignment of these domains resulting in finite bulk 
magnetization. The associated Tc is considerably lower 
than the one predicted by mean-field theory£3 This idea 
is also qualitatively depicted in Fig.|](c). In this scenario 
disorder enhances T*, but reduces Tc- 
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FIG. 3. MC time-averaged spin correlations per site vs. T 
(scaled by |S| 2 ) for the system shown in Fig.jj] for external 
fields h—0.0, 0.04, and 0.08. Correlations are measured either 
between sites i,j of the same cluster (filled symbols, cluster 
A or B), or between sites that do not belong to the same 
cluster. At high temperatures, the spin-spin correlations are 
approximately zero at all MC times. 

These percolation-like effects are also visible once a 
magnetic field is applied. As shown in the inset of Fig|l], 
small fields of (5-1 0T) are sufficient to introduce 

a robust change in the total magnetization of the sample 
at low temperatures. From Fig.[s| it can be concluded 
that this happens mostly through the mutual alignment 
of already preformed FM areas, whereas the magnetiza- 
tion of a given domain remains comparatively unaffected 
by those small fields. The susceptibilities should be even 
stronger for smaller values of x, especially in the insulat- 
ing regime. A pronounced magnetoresistance effect was 
indeed observed in lightly doped, insulating InAs and 
GaAs systems in a similar temperature regime.Ell The 
situation described here is reminiscent of the colossal ef- 
fects that can be observed in transition metal oxides. In 
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that context it was conjectured that disorder by chemical 
doping leads to cluster formation and concomitankJarge 
responses to small changes in various parameters. 113 
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FIG. 4. m vs. T for systems with different doping fractions 
p=0.05, 0.1, and 0.2, and a constant Mn density £=0.02 at 
zero external field. Weakened compensation leads to an en- 
hanced charge carrier rate, with extended FM domains. The 
inset shows the intercluster (A,B) spin correlations for p=0.20 
and p=0.05, demonstrating the influence of p on the correla- 
tions. 

The possible percolative nature of the magnetic tran- 
sition in DMS has been discussed recentlyE3, and it 
was proposed that the relevant analog is the problem 
of randomly placed spheres which overlap and, thus, 
form an infinite connected area once the radius of the 
spheres is larger than the percolation radius r per . In 
d=3, r per =1.33r , with r (ocrr -1 / 3 ) the interspin dis- 
tance (and percolation radius) in the periodic case. The 
MC data (r per «2-3 lattice units) agree well with this pre- 
diction. In this scenario an increase in x will lead to 
a smaller r per and, therefore, to an increase in Tc- At 



some point, however, 



will be smaller than the typical 



cluster size at T* , which makes T* the maximum criti- 
cal temperature possible. This temperature, although 
slightly lower, is in the range of the observed T™ ax and 
small changes in the parameters of Homs might be suf- 
ficient to result in higher transition temperatures. 

This can be demonstrated by altering the compensa- 
tion rate. Fig.^ shows m for the same impurity configura- 
tion as above, but with p=0.20 and p=0.05, respectively. 
As expected, a lowered compensation rate leads to en- 
hanced magnetism. The analysis of the spin-correlations 
(inset Fig.[|) shows that this is mostly due to an increased 
intercluster spin correlation, whereas the m vs. T curve 
for a given domain (not shown) is largely unchanged from 
the one at p=0.1. Thus, the increased number of charge 
carriers mainly leads to larger droplets, favoring the cre- 
ation of a coherent FM state. Additionally, the m vs. 



T curve acquires a more Brillouin-like look.y The oppo- 
site is true for the case of even stronger compensation, 
p=0.05, where a small reduction in Tc is observed. 

In summary, a model for DMS has been studied with 
MC methods at low temperatures. For the doping rates 
considered, the system can be described in terms of fer- 
romagnetic droplets. A change in parameters such as 
temperature or external fields can lead to the mutual 
alignment of these domains, while Tc is considerably re- 
duced compared to a mean-field analysis. The HF predic- 
tions of a disorder-enhanced transition point are not in 
agreement with nearly exact MC methods. The predicted 
Tc is somewhat lower than measured experimentally, but 
this value depends considerably on parameters that are 
not very well known. The inhomogeneous state is ex- 
pected to show large susceptibilities, e.g., a pronounced 
magnetoresistance effect. 
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